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^-\ • Flavour changing (FC) neutrino-matter interactions have been proposed as a 

lO . solution to the atmospheric neutrino anomaly. Here we perform the analysis of 

_i I the full set of the recent 52 kTy Super-Kamiokande atmospheric neutrino data, 

^ I including the zenith angle distribution of the contained events as well as the 

^} ' higher energy upward-going stopping and through-going muon events. Our 

V/-^ ' results show that the FC mechanism can describe the full data sample with 

^ . a Xmin — 44/33 d.o.f which is acceptable at the 90 % confidence level. The 

Q>^ i combined analysis confines the amount of FC to be either close to maximal 

G^ ; or to the level of about (10-50)%. 
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I. INTRODUCTION 

Neutrinos produced as decay products in hadronic showers from cosmic ray collisions 
with nuclei in the upper atmosphere ||l| have been observed by several detectors p|-§|,|[. 
Although the absolute fluxes of atmospheric neutrinos are largely uncertain, the expected 
ratio (/i/e) of the muon neutrino flux (z/^ + z/^) over the electron neutrino flux {ue + I'e) is 
robust, since it largely cancels out the uncertainties associated with the absolute flux. In fact, 
this ratio has been calculated [|I| with an uncertainty of less than 5% over energies varying 
from 0.1 GeV to 100 GeV. In this resides our confidence on the long-standing atmospheric 
neutrino anomaly. 

Although the first iron-calorimeter detectors in Frejus pf and NUSEX [^] reported a 
value of the double ratio, R(/i/e) = (/i/e)data/(At/e)MC5 consistent with one, all the water 
Cerenkov detectors Kamiokande [Q, 1MB [Q and Super-Kamiokande have measured 
R(/i/e) significantly smaller than one. Moreover, the Soudan-2 Collaboration, also using an 
iron-calorimeter, reported a small value of R(yu/e) 0], showing that the so-called atmospheric 
neutrino anomaly was not a feature of water Cerenkov detectors. 

Recent Super-Kamiokande high statistics observations indicate that the deficit in the 
total ratio R(/i/e) is due to the number of neutrinos arriving in the detector at large zenith 
angles. Although e-like events do not present any compelling evidence of a zenith-angle 
dependence, the /i-like event rates are substantially suppressed at large zenith angles. 

The z/^ — ^ Ut as well as the z/^ — >• z/^ P,pi| oscillation hypothesis provides a very good 
explanation for this smaller-than-expected ratio, which is also simple and well-motivated 
theoretically. This led the Super-Kamiokande Collaboration to conclude that their data 



provide good evidence for neutrino oscillations and neutrino masses |11]. However, alter- 
native explanations to the atmospheric neutrino data have been proposed in the literature 
including the possibility of neutrino decay |T2|, the violation of relativity principles p3|JT^ 



or the violation of CPT symmetry |T^. These explanations, however, have been challenged 
by the precise data of Super-Kamiokande on upward going muon events |]T6[ which allows to 



study the energy dependence of the neutrino survival (or disappearance) probability |jT^,|T8 
Based on such observations, both the possibility of an explanation of the anomaly in terms 
of neutrino decay [|1^ as well as the violation of relativity principles or the violation of CPT 
symmetry [|T^, have been disfavoured. 



In Ref. |^0[ an alternative explanation of the atmospheric neutrino data in terms of FC 
neutrino-matter interactions [^ was proposed, and it was shown that even if neutrinos 
have vanishing masses and/or the vacuum mixing angle is negligible, FC neutrino matter 
interactions could account for the Super-Kamiokande results on contained events providing 
an excellent description to the data, statistically as good as neutrino oscillations. The 



validity of this explanation was first questioned in Ref. ^j where the authors presented 
arguments against the FC neutrino-matter interaction solution on the basis of a fit to the 
up-going muons data from SuperKamiokande. 

In this paper we re-analyze the possibility of explaining the atmospheric neutrino 
anomaly by means of z/^ -^ Uj. conversion induced by flavour-changing neutrino-matter 
interaction which can be effective during the neutrino propagation in the Earth. We extend 
the analysis of Ref. |^ to the new set of Super-Kamiokande data by including also the 
up-going muon samples. 



II. MASSLESS NEUTRINO EVOLUTION WITH FC INTERACTION 

In our phenomenological approach we assume that the evolution equations which describe 
the Pfj, —>■ Ut transitions in matter may be written as 

' = v2gJ ^/rtnf'^^i, (1) 

where Va = ^a{r) (a = yU, r) are the probabihty amphtudes to find these neutrinos at a 
distance r from their creation position, \/2GFnf{r)ey is the z/^ + / -^ ^t + f forward 
scattering amphtude and \/2GFnf{r)el is the difference between the Ut — f and u^ — f 
elastic forward scattering amplitudes, with nf{r) being the number density of the fermions 
which induce such processes. 

The parameters e and e' contain the information about FC neutrino interactions. Such FC 




interactions may be accompanied by neutrino mass p2[ but this need not be the case p3| , p4 
One description would be to parametrize directly the FC interactions in terms of an effective 
four-fermion Hamiltonian. This could, for instance, arise by renormalization effects from 



the unification scale down to the electroweak scale in, say, supergravity models |2J]. An 
alternative more phenomenological way is to consider the existence of a tree-level FC process 
^a + f -^ '^f3 + f where / is an elementary fermion (charged lepton or quark). The interaction 
can be mediated by a scalar or vector boson of mass m and the neutrino-fermion coupling 
is generically denoted by Qaf (« is a fiavour index) and can be written as 

Since we are assuming vanishing neutrino masses, the anti-neutrino transitions z/^ —* Vr 
are governed by the same evolution matrix given in Eq. ([1|). For the sake of simplicity, 
we consider ep = e^ and ep = e'^,, which implies that we have only two free parameters in 
the analysis. Moreover, we set our normalization on these parameters by assuming that 
the relevant neutrino interaction in the Earth is only with down-type quarks. One could 
also assume that the incoming atmospheric neutrino has FC interactions off-electrons or 
equivalently, due to charge neutrality, off-up-type quarks. For simplicity, in the present 
analysis we consider only the case of interactions on down-type quarks. 

We have calculated the transition probabilities of z/^ — >■ v^ {p^ — > z^t) as a function of 
the zenith angle by numerically solving the evolution equation using the density distribution 

25| and a realistic chemical composition with proton/neutron ratio 0.497 in the mantle 



m 



and 0.468 in the core 26 



III. FITTING THE DATA TO THE FC HYPOTHESIS 

We have then used these probabilities to compute, as a function of the two parameters, 
Ey and e^, the theoretically expected numbers of events for the four sets of data reported by 
Super-Kamiokande: sub-GeV, multi-GeV, stopping muons and through-going muons. The 
expected number of contained events are computed by convoluting the probability with the 
corresponding neutrino fluxes (for which we use the Bartol calculations ||l|) and interaction 



cross sections and taking into account the experimental efficiencies as detailed in Ref. 0. 
For the up-going muon samples we obtain the effective muon ffuxes for both stopping and 
through-going muons by convoluting the probabilities with the corresponding muon fluxes 
produced by the neutrino interactions with the Earth. We include the muon energy loss 
during propagation both in the rock and in the detector according to p7| , p8| and we take 



into account also the effective detector area for both types of events, stopping and through- 
going. We compute the effective area using the simple geometrical picture given in Ref. 05 



Our final results show good agreement with the full MC simulation of the Super-Kamiokande 
collaboration in the Standard Model case (see the thick solid line in Fig. 3.) 

In our statistical analysis we adopted the technique PJ30|| of fitting separately the an- 
gular distributions of the fi- and e-like contained events (A^^ and N^,, i stands for sub-GeV 
and multi-GeV) and the up-going muon fluxes ($^, j = stopping, through-going). The 
expected number of events have been compared with the recent 52 kTy data reported by 
the Super-Kamiokande Collaboration [0] and the allowed regions in the {e^,, e^) plane have 
been determined from a x^ fit- In constructing the x^ function, we explicitly take into ac- 
count the correlation of errors, both of theoretical and experimental origin. Details on the 
definition of the correlation matrix for contained events can be found in Ref. , while the 
definition of the sources of errors and their correlations for the up-going muons fluxes are 



given in Ref. |31,30]. Here we simply summarize that we consider the overall normalization 
of the up-going muon fluxes to be affected by an uncertainty of 20% but in order to account 
for the uncertainties in the primary cosmic ray fiux spectrum we allow a 5% variation in 
the ratio between muon events in different energy samples. We further introduce a 10% 
theoretical error in the ratio of electron-type to muon-type events of the different samples. 
Other important source of theoretical uncertainty arises from the neutrino interaction cross 
section which at Super-Kamiokande ranges from 10-15 %. Uncertainties in the ratio be- 
tween different angular bins are treated, similarly to Ref. [0, by allowing a variation of 
5% times the difference between the mean bin cosines. With our definition we obtain, for 
instance, Xsm =122/(35 d.o.f) which means that the SM has a CL of 10~^^ ! Using this 
same x^ function for the case of oscillations we obtained allowed regions for masses and 
mixing angles very similar to those obtained by the Super-Kamiokande collaboration both 
for contained events as well as for upward going muons 0. 

In Fig. 1^ we show the contours of the regions allowed by the Super-Kamiokande data. 
The different panels of the figure refer to the fits performed over the different sets of data 
separately: (a) sub-GeV; (b) multi-GeV; (c) stopping muons; (d) through-going muons. The 
shaded areas are the regions allowed at 90% C.L., while the dashed and dotted contours refer 
to 95 and 99 % C.L., respectively. The condition used to determine the allowed regions is: 
X^ = xlnn + Ax^ where Ax^ = 4.6, 6.0, 9.2 for 90, 95 and 99 % C. L., respectively. 

The allowed regions for the contained events are, as expected, similar to the ones obtained 
in Ref. [^. The individual best fits now improve with respect to the analysis of the old data: 
Xl^^^ = 2.4/(8 d.o.f.) for the sub-GeV data (e^ = 0.196 and e'^ = 0.013) and to xLn = 6-4/(8 
d.o.f.) for the multi-GeV sample {e^ = 0.689 and e'^, = 0.284). The combination of the two 
sets of contained events leads to allowed regions which are analogous to the ones reported 



in the Ref. |]20| and which are not reproduced again here. The best fit point corresponds to 
e„ = 0.95 and e'^, = 0.084 with Xmin = 9.3/(18 d.o.f.). The goodness of the fit to the contained 
events in the FC-neutrino interaction scenario can be understood since the suppression 



of the expected event rates for contained events is the same for sub-GeV and Multi-GeV 
samples. The Super-Kamiokande collaboration has also measured the energy dependence 



of the up-down asymmetry for contained events [^ and this clearly indicates a strong 
energy dependence of the asymmetry for muon-like events in the momentum range 0.2 GeV 
< Pfj^ < 2 GeV. The asymmetry is consistent with zero at low momentum but significantly 
deviates from the expectation in the SM at higher momenta. One may naively expect that 
since the FC conversion mechanism is energy-independent it could be in contradiction with 
this measurement. However, one must notice that the average angle between the directions 
of the final-state lepton and the incoming neutrino ranges from 70° at 200 MeV to 20° at 
1.5 GeV, so that at low momenta the possible asymmetry of the neutrino flux is largely 
washed out. In Fig. ^ we plot, together with the Super-Kamiokande data, the momentum 
behaviour of the asymmetry in the FC-neutrino interaction scenario calculated for the best 
fit point to the contained event sample. As seen in the figure the agreement is excellent. 

In Fig. |T], we also show the regions which are allowed by the up-going muons samples of 
Super-Kamiokande. Panel (c) stands for stopping muons and panel (d) for the through-going 
sample. In the case of stopping muons, we see that, analogously to the contained events, the 
allowed region lies in the sector of the plane where the average survival probability is of the 
order of a half, which is what appears to be needed for explaining the data. Instead, in the 
case of through-going muons, the experimental data do not show such a strong reduction 
with respect to the theoretical calculations, and therefore the allowed region lies in the upper- 
left corner of the parameter space, which refers to a smaller transition probability. In both 
cases, the best fit point for each individual sample is good: Xmin = 1/(3 d.o.f.) for stopping 
muons {eu = 0.756 and e^ = 0.196) and Xmin = 10.3/(8 d.o.f.) for the through-going case 
(e^ = 0.081 and e'^ = 0.260). Both for the contained and for the up-going events, the best 
fits have the same level of statistical confidence as compared to the oscillation interpretation 
of the atmospheric neutrino data. This is shown in Table I, where we report the best fit 
values we obtain for the different data sets in the case of the FC-i/ interactions scenario and 



in the case of the neutrino oscillation scenario ||31|| . 

The allowed regions can be qualitatively understood in the approximation of constant 
matter density. The conversion probability in this case is 

Pi^.^'^r) = ^^^sin\^-VL), (3) 



where rj = \jAe1 + e'J'\/2GFnf. For Uf = Ud ^ ^rif, and e'^ < e^, the oscillation length in 
matter is given by 



Lose = — ^ 1.2 X 10=^ 
V 



2 mol/cc 



n„ 



^v 



km. (4) 



From Eq. (^ one can see that in order to have a relatively large transition probability, as 
required by the contained events and, also, by the stopping muons events, the FC parameters 
are required to be in the region e[, ^ ey and r] > t^ jR®- This last condition leads to a lower 
bound on ey. The island in Fig. |l|.(b) corresponds to ?7 ~ t^/ R®- 

The combination of the different data sets in a single x^-analysis is shown in Fig. |^. 
Panel (a) shows the combination of the full angular distribution of contained events with 



the total (unbinned) event rate of stop and through-going muons data, while panel (b) refers 
to the combination of all the angular distributions, including that of through-going muon 
events. In Fig. ^.(a) the information brought by the higher energy data is effective at the 
normalization level, since no information about their angular dependence is included. In this 
case the allowed region is still relatively large although the description is already worse than 
in the oscillation case as can be seen by comparing the corresponding xLin (20.4/(10 d.o.f.) 
for FC as compared to 9.6/(10 d.o.f.) for the oscillation scenario). This worsening is due to 
the fact that in the FC scenario the transition probability is energy independent while the 
data shows a smaller conversion for the higher energy through-going muon events. As seen 
in Fig. ^(b), when the angular information of both stopping and through-going muons is 
included in the data analysis, the description becomes even worse, mainly due to the angular 
distribution of the through-going data set. The allowed regions now form a set of isolated 
small 'islands'. The best fit point corresponds to e^, = 0.57 and e'^ = 0.45 acceptable at the 
90 % CL (xL„ = 44/(33 d.o.f.)). 

The behavior of the allowed regions can be understood by observing Fig. ^ where we show 
the angular distributions for the four cases: (a) sub-GeV; (b) multi-GeV; (c) stopping muons; 
(d) through-going muons. We show the distributions for the best fit point obtained from the 
combination of contained events with the total number of upward going muons Pi = {e^, e'^) 
= (0.17,0.28) and for the best fit point obtained from the analysis of the full data set P2 = 
(0.57,0.45). Although both points give a similar normalization to the up-going muon data 
samples, point Pi gives a better description to the angular dependence of the contained 
events, but it does not describe well the zenith angle distribution of the through-going muon 
events. As commented above, such point correspond to an effective FC-oscillation length of 
the order of the Earth radius. In this case we can see the imprints of the "oscillatory" sine 
behavior in the expected angular distribution of the up-going muon events. Such behavior, 
however, does not appear to be present in the Super-Kamiokande data, leading to a worse 
overall fit. In the case of multi-GeV contained events, this oscillatory behavior is averaged 
out due to the smaller angular resolution in the data and point Pi can give a good description 
of the data. On the other hand, point P2 gives a worse description of the contained events 
but fits better the shape of upward going muon data, with the exception of the last three 
angular bins of the through-going sample, where it does not produce a sufficient amount of 
through-going muons at angles < 6 < 20° below the horizon. 

IV. DISCUSSION 

What can we say about the required strength of the neutrino-matter interaction in or- 
der to obtain a good fit of the observed data? From our best fit results we obtain that 
the atmoshperic neutrino data can be explained if e^, > 0.4 and e'^, > 0.1. In terms of 
the neutrino-quark couplings introduced in Eq. (^ we see that for masses m ^ 200 GeV 
combinations of couplings grf ■ g^f and \grf\^ — IS'/^/P oi the order of 1 are needed. 

Model independent constraints on the e parameters can be extracted from their contri- 



bution to the Ufj^ neutral current cross section measured at low energy ||32|. These limits are 
stronger for interactions with quarks due to the better precision of the a'^^ ^ data as com- 
pared to the a'^Q g. We have estimated that these data constrain e ^ (9(0.1)-C(1), depending 
on the fermion / coupled to the neutrino. No limit on the e' parameter can be obtained from 



these measurements. From this, we see that the FC mechanism is somewhat disfavoured, 
although not strictly ruled out. Additional limits can be obtained from the non-observation 
of lepton flavour violation in r decays. Indeed, if SUl{2) symmetry is assumed 0] these 
limits are rather stringent. However, they strongly depend on the amount of SUi{2) vio- 
lation present in the model. For the purpose of illustrating this explicitly, let us consider 
a supersymmetric model with broken i?-parity as a way to parameterize the FC neutrino- 



matter interaction [P3[. In this case the FC z/^-matter interactions are mediated by a scalar 
down-type quarks, dj, so that we need only to check the couplings where a c?-quark and a 
11- or r-neutrino is involved, ie gid ~ A^^-^, i = 2,3. The X'^^^ are the coupling constants in 
the broken i?-parity superpotential X'-ji^LiQjDl, where L, Q and D are standard superfields, 

and 4:\/2Gf£u = iXz-^sji-^sji/^jl- "^^^ rnost stringent limit to the values of the relevant FC 

j 
quantities comes from limits on the FC tau decay BR(r^ — > p° + n^) < 6.3 x 10^^ which 

implies that lX!'^3ii'^2ji(100 GeV/muj)^| < 3.1 x 10^^. This constraint can be satisfied 

j 
by cancellation between the contributions from the up-squarks exchange of the third and 

second generation, although a certain degree of fine-tuning is needed between their masses 
and couplings. Notice that the e,y may still be large as long as the cancellation does not 
occur for the down-quarks. This can be achieved by a small splitting between the up and 
down-squark masses without conflicting with the limits from Ap. 



V. CONCLUSIONS 

In summary, in this paper we have re-analyzed the possibility of explaining the atmo- 
spheric neutrino anomaly by means of i/^ -^ v^ conversion induced by flavour-changing 
neutrino-matter interaction which can be effective during the neutrino propagation in the 
Earth. We extend the analysis of Ref. [^ to the new set of Super-Kamiokande data by 
including also the up-going muon samples. Our results show that that flavour changing 
i/^-matter interactions are able to describe the full set of data of Super-Kamiokande on at- 
mospheric neutrinos at the 90 % CL. The agreement between the data and the calculated 
events for the Super-Kamiokande detector is notably good for the individual sets of data 
collected by Super-Kamiokande, with a confidence level as good as for the oscillation hy- 
pothesis. When the data are combined together, in particular once the upward-going muon 
zenith-angle distribution is included in the analysis, the z/^ oscillation provides a much bet- 
ter description. However, the FC mechanism is compatible with the data at the quoted 90% 
level of statistical confidence. The worsening of the fit which occurs when the through-going 
muons sample is included is partly due to the fact that in the FC scenario the transition 
probability is energy independent while the data shows a smaller conversion for the higher 
energy up-through-going muon events. The ensuing result is that the expectations from 
FC-neutrino interaction for neutrinos arriving mainly at angles above 20 degrees below the 
horizon do not reproduce the experimental data. 

The amount of FC neutrino interactions required by our combined analysis in order to 
fit the data is somewhat large, either close to maximal or at the level of about (10-50)%. 
Although significant, FC at this level in the neutrino sector can be accomodated in speficic 
models without conflicting with existing experimental limits. These realizations of large 



neutrino FC may require some degree of fine tuning, like for instance in broken R-parity 
super symmetric models, but are theoretically viable. 

The above FC mechanism can be also tested at future Long Baseline experiments. From 
Eq. @, using ng ~ 2 mol/cc, we can predict that for e ~ e' ~ 1 (0.1) the planned K2K 
experiment [^] should obtain P^i'^ -^ Vt) ~ 0.35 (0.004) while for MINOS |^ one finds 
P{vt, -^ Vr) ~ 0.75 (0.04). 

In conclusion, we have to comment that although FC interactions would be eventually 
ruled out as the only source of the modification of the atmospheric neutrino predictions with 
respect to those of the Standard Model they could still be there at some level, even if the 
data would admit a very good interpretation in terms of standard v^ —>■ Vj. oscillations. This 
is theoretically not an ad hoc assumption, since in many theoretical models neutrino masses 
naturally co-exist with FC-neutrino interactions. 

On the other hand, one may turn the argument the other way around: should the 
atmospheric neutrino anomaly be explained in terms of neutrino oscillations, then it will be 
possible to use the non-observation of an additional effect in the atmospheric neutrino data 
in order to impose new model independent limits on the strength of the e and e' parameters, 
as can be foreseen by looking at Fig. |l](d). 
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TABLES 



Data 


d.o.f 


2 
XminFC 


v2 
''^minOsc 


sub-GeV 


8 


2.4 


2.4 


multi-GeV 


8 


6.4 


6.3 


contained 


18 


9.3 


8.8 


stopping-^ 


3 


1. 


1.3 


through-going-^ 


8 


10.3 


10.4 


contained -|- total up-^ 


10 


20.4 


9.6 


contained -|- angular up-fi 


33 


44. 


23.5 



TABLE L Xmin obtained for several data combinations in the framework of FC-i^ interactions 



as compared to the case of the neutrino vacuum-oscillation scenario. 
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FIG. 1. Allowed regions for e,y and e'^ in the FC massless-neutrino scenario for the different Su- 
per-Kamiokande data sets: (a) sub-GeV, (b) multi-GeV, (c) stopping muons and (d) through-going 
muons. The best fit points for each case are indicated by stars. The shaded area refers to the 90% 
C.L.. while the contours stand for 95% and 99% C.L. 
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FIG. 2. Expected up-down asymmetry in Super-Kamiokande for the FC-neutrino scenario 
as a function of the muon momentum for fully contained ^-like events, compared to the Su- 
per-Kamiokande experimental data. 
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FIG. 3. Allowed regions for €i, and e^ for the combination of the Super-Kamiokande data sets: 
(a) the binned contained events are combined with total (unbinned) up-going events; (b) binned 
contained and up-going events. The best fit points for each case are indicated by stars. The shaded 
area refers to the 90% C.L. while the contours stand for 95% and 99% C.L. 
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FIG. 4. Best-fit zenith angle distributions in the massless- neutrino FC scenario. The thick-sohd 
lines correspond to the calculation in absence of new physics. The dotted lines correspond to the 
best fit point obtained by the analysis of the contained events combined with total (unbinned) 
up-going events. The thin-solid line is for the best fit point of the combined analysis of contained 
and up-going muon events. The 52 kTy Super-Kamiokande data are indicated by crosses. 
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